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^ ' We have investigated plasma-surface interactions with molecular dynamics (MD) simulations.^'^ It, 

however, is high cost computation and is limited to simulations for materials of nanometer order. 
In order to overcome the limitation, a complementary model based on binary collision approximation 
(BCA) can be established. We employed a BCA-based simulation code ACAT^ and extended to handle 

C^ ' any structure involving crystalline and amorphous. The extended code, named "ACVT", stores all 

positions of projectile and target atoms and velocities of recoil atoms, so it can be combined with the 
'^ ■ MD code. It also holds the potential to reproduce channeling phenomena. Thus it is expected to be 

Q ' useful for evaluation of channeling effects. 

KEYWORDS: binary collision approximation, plasma-surface interactions, sputtering, graphite, hy- 
►^ . drogen 

m 

. ■ 1. Introduction 

in 

<^ ' In nuclear fusion devices plasmas of hydrogen isotopes exist and contact material surfaces. 

o : 

The "divertor configuration" is employed in order to control impurities produced by impacts 
of plasmas to the surface and to reduce the heat load to plasma facing materials. In the 
^ ■ configuration divertor plates, whose potential constituent includes carbon and tungsten, are 

C^ ' installed. Understanding of the divertor physics and designing an appropriate configuration is 

essential for the establishment of nuclear fusion reactor. These require knowledge on plasma- 
surface interactions (PSI). Thus we have investigated interactions between hydrogen atoms 
and carbon material, such as graphite, with molecular dynamics (MD) simulations.^'^ 

The molecular dynamics simulation code solves equations of motion for all particles un- 
der the modified Brenner reactive empirical bond order (REBO) potential.^ Evaluation of 
the REBO potential requires consideration of effects from multiple particles, and thus it is 
generally high cost computation. This limits the applicable material scale length of the MD 
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simulation to about order of nanometer and the energy range to about order of kilo electron 
volt. 

In order to investigate PSI with numerical simulation we have to overcome the limitation 
contained by MD simulation. A complementary model can be established based on binary 
collision approximation (BCA). There exist intensive PSI related works with binary collision 
based monte carlo simulations.^"^ The binary collision approximation simplifies interactions 
between material elements and reduces them to the sequence of the binary collisions. A benefit 
of the model is that it is rather simple and requires less computing resources than MD model. 
It, however, holds a limit of application on lower energy region. 

We think that a hybrid simulation of molecular dynamics and binary collision approxima- 
tion is promising. In the integrated simulation code the part of binary collision approximation 
covers higher energy region, and the part of molecular dynamics governs lower energy region 
and solves only the vicinity of the projectile and recoil atoms. The threshold energy between 
the two regions should be determined according to a validity condition of the binary collision 
approximation, which is roughly estimated as 200 eV. 

In this paper we extend an existing binary-collision-approximation-based simulation code 
ACAT^ in order to combine MD simulation code and BCA based one. 

2. ACAT code and its extension code ACVT 

We employ a binary-collision-approximation (BCA) based simulation code ACAT. The 
ACAT code was developed to simulate atomic collisions in an amorphous target within the 
framework of the binary collision approximation. Projectile particles are traced through binary 
collisions. Target particle, with which projectile collides, is randomly distributed in each unit 
cell whose size Rq = N^^'^, where N is the number density of the target material. In terms 
of randomly distributed target particle, this code employs the monte carlo method and aims 
for the atomic collisions in amorphous target. 

Figure 1 depicts the trajectory of two particles interacting according to a conservative 
central repulsive force. The scattering angle in the center-of-mass system (CM-system) is 

f°° 1 
e = 7T-2b -^r-^dr, (1) 

where, 



62 V{r} 



5(0 = ^/1-12-^' (2) 



b is the impact parameter, Er = EQmi/{mi + 1712) is the relative kinetic energy, Eq is the 
incident kinetic energy of the projectile, V{r) is the interatomic potential, rg is the solution 
of g{r) = 0, nil and 1112 are the mass of the projectile and the target atom, respectively. 

The trajectories of particles are approximated as the asymptotes of them in the laboratory 
system (L-system). So they consists of linkage of straight-line segments. The starting point 
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of the projectile and the recoil atom after a collision is given by Axi and Ax2, which are the 
shifts from the initial position of the target atom shown in Fig.l: 

_ 2r + (^-l)6tan(e/2) 



Ax2 = 6tan(e/2)- Axi, (4) 

where 



and the mass ratio A = 1112/1711. 

As the interatomic potential V{r), the Moliere approximation to the Thomas-Fermi po- 
tential^ is employed: 

V{r) = ^^^^$(r/a), (6) 

$(x) = O.Sbe-^-^'-" + 0.55e-i-2^ + 0.10e-^•°^ (7) 

where a is the screening length, and Zi and Z2 are the atomic numbers of the projectile and 
the target atom, respectively. 

The procedure of searching the collision partner is as follows: We define a unit vector Cp 
in the direction of a moving projectile and a step length Ax. Here, the notation "projectile" 
means not only the incident particle but any recoil atom. The position of the projectile in 
n-th step , which is initially R^^', is i?'"-* = R^^' + nAxCp. If the unit cell involving R^^' is 
different from that of the initial position R^ ' , a target atom is produced in the new unit cell 
by use of four random numbers. Three random numbers are for location and one is for the 
kind of target atom. This target atom is the partner in the next collision. 

The impact parameter b is given as 

b=\{RA-R)xep\, (8) 

where Ra is the position of the target atom. Figure 2 depicts the situation described above. 
We can obtain the trajectory of particles by the eqs. 1-8. 

Three parameters are defined: bulk binding energy E^, displacement energy E^^, and 
minimum energy E^.. Let us consider a collision from which the original projectile emerges 
with kinetic energy Ei after transferring kinetic energy T to the target (collision pair) atom. 
The target atom is displaced when its kinetic energy exceeds the displacement energy E^^. If 
T > E'd is satisfied, the target atom is added to the collision cascade with the kinetic energy 
E2 = T — Eb- Each projectile is traced while its kinetic energy exceeds a given minimum 
energy Ec- 

We have extended the ACAT code described above as follows: While a target atom is 
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produced randomly within a cubic collision cell in the original code, our extended version 
employs rectangular parallelepiped cells and positions of target atoms are given as initial 
condition. The positions of projectile and target atoms are stored in computer memory, which 
enables us to treat structured target materials as well as amorphous ones. In the present 
version it is assumed that each cell contains zero or one target atom. Thus we can obtain 
the position of target (collision pair) atom from the stored data when the collision cell is 
specified. We call the extended version "ACVT" code, which stands for atomic collisions in 
any structured target. The notation V implies that the code can handle any structure involving 
crystalline and amorphous. 

If the kinetic energy of each projectile T becomes less than a predetermined minimum 
energy Ec, for example 200 eV, tracing of the projectile is stopped and its location and 
velocity are stored, which can be used for initial condition of a succeeding molecular dynamics 
simulation. 

3. Application of ACVT code to hydrogen injection into graphite 

A graphite with the size of 32lA (W) x 347A (D) x 335A (H) is setup as a target material, 
which consists of 100 layer graphenes. This material is put into a simulation domain, which is 
divided into cells whose size is (v3?'o/2) x (3ro/2) x (/iq), as shown in Fig. 3. Here rg = 1.4463A 
is the distance between neighboring carbon atoms and /iq = 3.3480A is the distance between 
graphene layers. Then each cell contains one carbon atom and the simulation domain includes 
256 X 160 X 100 cells. 

As an example problem, we inject a hydrogen atom with 1.0 keV into the target graphite. 
Incident angle is fixed to degree, that is, normal to the target surface. Incident point is 
varied for each trial. 

Figures 4(a) and 4(b) show trajectories of two cases of an example problem. Differences 
of input parameter between two cases are only on the hydrogen incident point. Black solid 
circles and gray ones denotes initial points and final points, respectively. 

The top black point corresponds to the hydrogen incident point. The hydrogen atom 
(projectile) moves along the curve with small angle scatterings. In a small angle scattering, 
energy transferred from the hydrogen atom to a target carbon atom is so low that the target 
cannot be displaced. When the transferred energy exceeds the displacement energy E'^, a 
recoil carbon atom, which initially exists at a black solid circle in the figures, is also traced. 
If the kinetic energy of an atom becomes less than the given minimum energy Ec, tracing the 
particle is stopped. The final position is depicted as gray solid circles in the figures. 

In figures 4(a) and 4(b), we can see some clusters where a few solid circles gather. In 
each cluster, a recoil carbon atom collides with another carbon atom, transfers its own kinetic 
energy to the collision pair, and the collision pair recoils. That is, collision cascades occur in 
each cluster. 



4/10 



Jpn. J. Appl. Phys. Regular Paper 

In figure 4(b), the curve of trajectory in upper lialf is almost straight. In this case the 
projectile of hydrogen moves straight about 75A along the interstitial, which may be regarded 
as a kind of channeling. This kind of trajectory is not realized in the original ACAT code 
because it assumes amorphous target. In figure 4(a), the projectile trajectory just before 
its stopping lies horizontally with fluctuations. In this region the projectile travels between 
graphene layers. This kind of trajectory is also newly reproduced by our extension of the 
ACAT code. We expect that the ACVT code is useful for evaluation of channeling effects. 

4. Conclusion 

We have extended the ACAT code, which stands for atomic collisions in an amorphous 
target and is based on a binary collision approximation, to handle any structure involving 
crystalline and amorphous. It is shown that new types of trajectories can be reproduced by 
use of the extended version, named "ACVT" code. The name "ACVT" stands for atomic 
collisions in any structured target, and the notation V implies that the code can handle any 
structure involving crystalline and amorphous. 

Now we develops a hybrid code of our molecular dynamics (MD) simulation code and 
binary-collision-approximation (BCA) based code "ACVT" using a multiple-program multiple- 
data (MPMD) approach. ^"^ As the BCA-based code is relatively lightweight, the resultant 
hybrid simulation code is expected to simulate carbon materials of submicron order. 

The ACVT code itself holds new features and is promising for evaluation of channeling 
effects. 
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Fig. 1 . The trajectory of two particles interacting according to a conservative central repulsive force in 
the laboratory system (L-system). The positions of the projectile and the target atom correspond 
to the apsis of the collision.^ 
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Fig. 2. The procedure of searching the collision partner:^ A unit vector Cp is in the direction of a 
moving projectile. Projectile starts from i?^ ' , moves straight in Cp-direction with the step length 
Ax. If the unit cell involving the projectile after n-th step i?^"-* is different from that of the initial 
position R" ' , a target atom Ra is randomly produced in the new unit cell. The variable b denotes 
the impact parameter. 
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Fig. 3. The configuration of tlie target materials of graphite and the construction of simulation cell 
in (a) cross-sectional view, (b) vertical cross-sectional view. Dots and their linkage denote carbon 
atoms and bonds, rg = 1.4463A and 3.3480A correspond to the length of the bond and the distance 
between graphene layers, respectively. In fig. (a), the bonds in black and ones in gray belong to 
another graphene layer. Rectangles denote cells for ACVT simulation. 
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Fig. 4. Trajectory of hydrogen injected into graphite and trajectory of recoil carbon atoms. Injection 
energy of hydrogen is 1.0 keV and incident angle is degree. Incident point differs in (a) and (b). 
The size of boxes are (a) lOOA x 60A x 140A and (b) 250A x 120Ax 200A. Black sohd circles and 
gray ones denote initial points and final points of the injected hydrogen atom and recoil carbon 
atoms. 
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